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Dynamic correlations of quantum magnets provide valuable information. But especially at fi¬ 
nite temperature, many observations, e.g., by inelastic neutron scattering experiments, are not yet 
quantitatively understood. Generically, the elementary excitations of gapped quantum magnets 
are hardcore bosons because no two of them can occupy the same site. The previously introduced 
diagrammatic Bruckner approach dealt with hardcore bosons of a single flavor at finite tempera¬ 
ture. Here, this approach is extended to hardcore bosons of several kinds (flavors). The approach 
based on ladder diagrams is gauged with so far unknown rigorous results for the thermal occupation 
function of multi-flavored hardcore bosons in one dimension with nearest-neighbor hopping. For 
low temperatures, it works very well. Furthermore, we study to which extent the approach is a 
conserving approximation. Empirical evidence shows that this is true only in the single-flavor case. 

PACS numbers: 75.40.Gb, 75.10.Pq, 05.30.Jp, 78.70.Nx 


I. Introduction 


A. General Context 

Understanding the finite temperature effects of correlated 
quantum systems is a complicated problem. Even at zero 
temperature the complex interplay of different interac¬ 
tions can lead to a wide range of interesting phenom¬ 
ena, such as bound states, fractionalization, quantum 
phase transitions (QPTs) or Bose-Einstein-Condensation 
|l|. Even more, at finite temperature, interesting fea¬ 
tures emerge from the interplay of thermal and quantum 
fluctuations. Especially dynamic correlations are directly 
affected by microscopic interactions as well as by thermal 
effects. 

From a theoretical point of view, it is often advanta¬ 
geous to describe the low-energy physics of quantum sys¬ 
tems by an effective model which includes the motion 
and interactions of quasi-particles of which the number 
is conserved or almost conserved. There are numerous 
ways to derive effective models, an incomplete list com¬ 
prises renormalization approaches pi Q, continuous uni¬ 
tary transformations (CUTs) [41-1 1 Oh and variational ap¬ 
proaches, see for instance Refs, EH [l2j for variational 
approaches based essentially on intuition and Refs. El- 
1151 for variational approaches based on matrix product 
states (MPS). 

Variational approaches and CUTs generically provide ef¬ 
fective models at zero temperature. Thus they do not 
address the thermal effects directly which renders quan¬ 
titative comparisons with experiments difficult. But if 
the effective model is formulated in second quantization 
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it can be evaluated also at finite temperatures by addi¬ 
tional theoretical tools. Such tools are for instance ex¬ 
act diagonalization, thermodynamic Bethe ansatz, or di¬ 
agrammatic approaches. 

Hardcore bosons are one generic type of elementary exci¬ 
tations in quantum lattice systems [§). They behave like 
normal bosons on different lattice sites, but at one site 
only a single excitation is allowed at maximum. In con¬ 
trast to free bosons or free fermions, the thermal proper¬ 
ties of hardcore bosons are unknown and not trivial even 
if they are not subject to further interactions. 

The motivation to study hardcore bosonic excitations 
originates from quantum magnets HQ- A variety of 
elementary hardcore bosonic excitations exist in quan¬ 
tum magnets, e.g., spin-flip excitations in polarized sys¬ 
tems, domain walls in one-dimensional systems with de¬ 
generate ground states, or triplon excitations in strongly 
dimerized systems. Experimentally, inelastic neutron 
scattering (INS) provides a convenient tool do study the 
effects of finite temperatures on dynamic correlations 
Pi , namely by probing the spin dynamic structure factor 
(DSF) 

I 

S a (u,Q) = - Y, eiUt e- m - l) (S?(t)S?) (1) 
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where L is the number of sites. 

In this article we focus on gapped hardcore boson sys¬ 
tems. Quantum antiferromagnets built from dimers of 
5 = 1/2 spins on the lattice are of special interest. Exam¬ 
ples include spin chains 7|. ITol Ip . spin ladders Esm- 
l24j . coupled ladders gml and dimers coupled in sev¬ 
eral spatial directions MMil. This underlines that 
the topic of the present article is of current interest and 
relevant to experiments although it focuses on conceptual 
issues and benchmarks of the approach. 

For gapped quantum magnets the generic DSF spectrum 
consists of a dominant single quasi-particle peak and 
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weaker multi-particle continua. At finite temperature, 
the single-particle peak is broadened due to thermal scat¬ 
tering. In recent years, the energy resolution of INS ex¬ 
periments improved significantly. This makes it possible 
to investigate the position and the shape of the response 
as function of temperature. 

Various theoretical approaches have been applied to de¬ 
scribe the effects of finite temperature on the dynamics 
in quantum magnets. One of them is complete diag- 
onalization to compute the Lehmann representation of 
the DSF |33j, [HI- For one-dimensional systems, Essler 
and co-workers developed a variety of approaches based 
on integrabilit y, B ethe-Ansatz equations and fermioniza- 
tion in Refs. [35M38I ] . Certain thermal effects such as 
the narrowing of the dispersion and shifts of spectral 
weight can already be captured by mean-field approaches 
fjjl . The combination of integrability and density-matrix 
renormalization is a very useful strategy where appli¬ 
cable [39], An advanced numerical approach is based 
on density-matrix renormalization formulated in matrix 
product states in the frequency space [13 ■ 

The above listed methods rely essentially on the one- 
dimensionality of the system under study. In principle, 
the complete diagonalization can be applied to any di¬ 
mension. But the tractable system size measured in the 
linear extension becomes very small for dimensions higher 
than one. Thus, recent finite temperature experiments 
on 3D materials [30| can not be described by these ap¬ 
proaches. A 1 j z expansion with z being the coordination 
number by Jensen [311 ] agrees remarkably with the exper¬ 
imental data. The underlying idea relies on a description 
of the scattering of the elementary excitations with an 
effective medium 13T1 PiTil . The caveat of the approach is 
that it is unclear how it can be systematically extended 
beyond the linear order in \/z. Already in this linear 
order, non-physical double poles have to be eliminated 
by a plausible assumption in order to preserve causal¬ 
ity. 


B. Diagrammatic Approach 

In 2014, we introduced a diagrammatic perturbative ap¬ 
proach based on a low-temperature expansion which does 
not require a particular dimensionality [43 ]. The system¬ 
atic control parameter of the approach is the low density 
of thermally excited hardcore bosons which is propor¬ 
tional to exp(—/3A) in a gapped system where /? is inverse 
temperature and A the energy gap. 

The method works conceptually in any dimension; it di¬ 
rectly addresses frequency and momentum space in the 
thermodynamic limit. Therefore finite size effects and 
the ill-defined continuation from imaginary Matsubara 
frequencies to real frequencies are avoided. Furthermore 
the approach lends itself for effective models at T = 0 
derived by other methods. 

A key ingredient for the applicability of the advocated ap¬ 
proach is that the system conserves the number of quasi¬ 


particles (elementary excitations), i.e., there are no finite 
anomalous Green functions and the vacuum is the ground 
state of the system. 

One may think that this property contradicts the exis¬ 
tence of quantum fluctuations in the original model. But 
this is not the case. The idea is to treat the quantum fluc¬ 
tuations beforehand by another approach, based on vari¬ 
ation or on unitary transformations. The effective model 
obtained should conserve the number of quasi-particles 
with the ground state being their vacuum. On top of this 
effective Hamiltonian, we apply the Bruckner approach 
to treat the thermal effects in the spectrum. 

The method is motivated by the since long known Briick- 
ner theory in nuclear matter [33 ] and is discussed as 
the T-matrix approximation in Ref. [46j. But it has also 
been employed for spin systems at zero temperature by 
Sushkov and co-workers [47], !48|. In these applications, 
however, the quantum fluctuations were treated diagram- 
matically since the Bruckner idea was not applied to a 
particle-conserving model [33 ] which makes its applica¬ 
tion non trivial even at zero temperature. 

The Bruckner approach takes into account that the ele¬ 
mentary excitations in quantum paramagnets are hard¬ 
core bosons, i.e., a site can be occupied at maximum by 
a single excitation. This is enforced by an infinite on-site 
repulsion. Note that at zero temperature the hardcore 
constraint is irrelevant for the single-particle dynamics if 
the Hamiltonian is particle conserving. At finite temper¬ 
ature, however, the hardcore scattering provides the pri¬ 
mary source for the loss of coherence of the elementary 
excitations. Then, the single particle is scattered from 
the thermally excited background of particles. 


In a first article [44j] we introduced the Briickner ap¬ 
proach at finite temperatures and tested it for hardcore 
bosons with a single flavor in one dimension. The advan¬ 
tage of this model was that it can be solved using the 
Jordan-Wigner transformation [5(|. It turned out that 
the Bruckner approach based on ladder diagrams works 
very well. In particular, we empirically found that the in¬ 
tricate sum rules of hardcore bosons were conserved even 
at high temperatures at and beyond the energy gap A 
within numerical accuracy. Furthermore, the approach 
yielded a very good approximation for the exact thermal 
occupation for all temperatures. This has been recently 
confirmed and quantified even down to vanishing gaps 


51]. The asymmetric line-shape broadening and band¬ 


width narrowing found in other theoretical studies EH- 
l38j and in experimental studies, e.g., in Ref. [2(j, are 
retrieved. 


In the present paper, we extend the finite-temperature 
Bruckner approach to systems with multiple flavors, i.e., 
there are hardcore bosons of different kinds. But each 
single site can still only be occupied by a single boson. 
This situation is relevant wherever the local site is repre¬ 
sented by a more complex quantum system with a Hilbert 
space dimension larger than two. A generic example are 
systems built from coupled dimers, each consisting of two 
spins S = 1/2. In these systems the elementary excita- 
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tions are mobile S = 1 triplets, which we call triplons 
to distinguish them from magnons in ordered magnets 
[hd ]. For full SU(2) symmetry the triplons are threefold 
degenerate, i.e., they are hardcore bosons with three fla¬ 
vors. 

The goal of our study is to investigate the multi-flavor 
extension of the Briickner approach and to gauge the 
method against exact results. Note that a benchmarked 
method to calculate finite temperature correlations ap¬ 
plicable in any dimension is essential to describe and to 
explain experimental observations in quantum magnets. 
We study the performance of the method with a special 
focus on sum rules and discuss to which extent our ap¬ 
proach represents a conserving approximation for hard¬ 
core particles. 

We stay here with the one-dimensional case because it 
allows us to compare the diagrammatic approach to rig¬ 
orous results. But it must be emphasized that there is 
no conceptual obstacle for the application of the diagram¬ 
matic approach in higher dimensions as well. 

In order to have rigorous results at our disposal we de¬ 
rive the thermal occupation function for hardcore bosons 
in one dimension exactly for a nearest neighbor hop¬ 
ping- 

The paper is organized as follows: In Section Hi Al wc in¬ 
troduce the multi-flavor hardcore boson model and dis¬ 
cuss some of its properties. In the following Section 
III Bl we derive the exact thermal occupation for a special 
case of hardcore bosons in one dimension with nearest- 
neighbor hopping. In Section Hill we extend the Briickner 
approach on the basis of ladder diagrams to the multi¬ 
flavor case. Next, in Section HVl we evaluate the results 
of the diagrammatic Bruckner approach on the basis of 
the findings of the preceding sections. In Section [V] we 
perform various calculations at finite repulsion, i.e., the 
bosons repel each other only with a finite energy cost, 
to clarify to which extent conserving approximations ex¬ 
ist for hardcore bosons. Finally we conclude in Section 

ED 


{ i , d, di, d, 2 , ^ 3 } denote site indices which need not be re¬ 
stricted to one dimension, {a, (fi, 7 , £} denote flavor in¬ 
dices, are hopping matrix elements, are in¬ 

teraction vertices and b\ a ,b ia denote hardcore bosonic 
creation and annihilation operators. The hopping ma¬ 
trix elements must fulfill the hardcore constraint, i.e., 
d 0 , di 0 , d 2 7 ^ d 3 . 

The model @ does not include terms that change the 
number of particles. These terms would induce addi¬ 
tional quantum fluctuations and the vacuum would not 
be the ground state of the system. It is understood that 
particle number violating terms have already been ac¬ 
counted for by other methods, e.g., based on unitary 
transformations or on variational principles. Hence we 
do not neglect quantum fluctuations, but the Hamilto¬ 
nian m describes the hopping and the interactions of 
conserved effective particles. We aim at the systematic 
description of thermal effects in spectral properties for 
the above effective Hamiltonian. 

In Ref. 0 we considered the single-flavor case of a chain 
with nearest-neighbor hopping without any additional in¬ 
teractions. Here, we address the consequences of intro¬ 
ducing several flavors. 

The hardcore operators in second quantization fulfill the 
algebra, 
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Transforming the hardcore commutator into Fourier 
space yields 


bk,a,b\, ^ ( 1 b\, +gsi bk+q, 1 


b k'+q,<t> bk +<h a 


( 4 ) 


II. Model and Exact Results 
A. General Properties 

In this section, we introduce the multi-flavor hardcore 
boson model and discuss some of its properties such as 
commutators, spectral functions, and sum rules. 

The general Hamiltonian of conserved hardcore bosons 
reads 

H 0 = E 0 + J2 lA+d,a + hx -) 

i,d,ot 

^ V d 1 ,d2,d 3 U i,a u i+d 1 ,</> U i+d 2 ,i U i+d 3 ,£T ■■■’ \ A > 

i, di, d- 2 , d 3 
a. </>, 7> £ 

where the dots stand for possible higher particle-number 
interactions. Here Eq denotes the ground state energy, 


The quantity of interest is the single-particle spec¬ 
tral function depending on momentum p and frequency 
w 


A a [p,u) = — lim 

7T iu) v — >-a;+iO+ 


p 

Im J d re^ T ^= £ e~™G a {j, r), (5) 


where G a (j, r) is the single-particle temperature Green 
function 

G a {j, t) = - (T{4(-ir)6 0iQ (0)}) . (6) 


We use the Matsubara frequencies = 2vi t//3 and the 
inverse temperature /3 = 1/T. For hardcore bosons the 
spectral function is a non-trivial quantity because it can 
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not be calculated exactly, even in the absence of any addi¬ 
tional interactions beyond the hardcore constraint. The 
hardcore constraint itself represents a strong interaction 
inducing intricate correlations. 

The spectral function is related to the dynamic structure 
factor (DSF) by means of the fluctuation-dissipation the¬ 
orem 

■S a (P, “)= 1 _ \- Pu [^°(P. tt) + A a (p, -u )]. (7) 

The DSF is directly accessible by scattering experi¬ 
ments. 

Sum rules, i.e., the integrals of the spectral function over 
frequency, often provide a convenient way to evaluate the 
accuracy of various approximations. For normal bosons, 
the integration of the spectral function over frequency 
yields a temperature independent constant of unity. But 
for hardcore bosons this is no longer true. The general 
sum rule states that the equality 

OO 

J A a (p, w)dw = ( [b p>a , bl a ]) , (8) 


holds. Evaluating Eq. m with the commutator in 
Fourier space leads to 

OO 

J A a (p,uj)du) = 1 — (1 + N{)n(T), (9) 

— OO 

where Nf denotes the number of different flavors, e.g., 
iVf = 3 in the triplon case. The thermal occupation 
is a function of temperature and 
independent of the flavor if the different flavors are de¬ 
generate. 

At zero temperature, the thermal occupation vanishes 
and the standard boson sum rule is recovered. At infi¬ 
nite temperature, the thermal occupation function yields 
n(T ) = 1/(1 + Nf), so that the integration of the spec¬ 
tral function over frequency yields zero. In between, the 
thermal occupation is a non-trivial quantity. 


B. Exact Thermal Occupation in One Dimension 

In this subsection, we calculate the thermal occupation 
of multi-flavored hardcore bosons with nearest-neighbor 
hopping along a chain. We consider a chain of L sites of 
hardcore bosons with a nearest neighbor hopping Hamil¬ 
tonian without any additional interactions V 


The energy gap is given by A > 0 while W > 0 is the 
band width of the dispersion 

w{k) = A + — [1 — cos(fc)] . (11) 


Periodic boundary conditions are implicit and k is the 
momentum of the excitations. The minimum A of the 
dispersion is found at k = 0. The ground state of the 
system is given by the vacuum state. 

Since we only allow nearest neighbor hopping, the bosons 
cannot skip over each other. This implies that a given 
sequence of the flavors along the chain is conserved in the 
dynamics induced by H in (fTUl) . Thus the flavor degree 
of freedom and the hopping dynamics can be considered 
to be completely decoupled for a fixed number of bosons. 
The flavor degree adds a certain degree of entropy for 
given number N of bosons. The hopping of the hardcore 
bosons acts as if all bosons have the same color. 

To calculate the thermodynamics of the flavorless hard¬ 
core bosons we apply a Jordan-Wigner transforma¬ 
tion 


Cj = exp(? Ti^blb^bj 
i<j 

bj = exp(-7 TiJ2 c Ui)cj 

i<j 


c) = exp(— to iJ2blbi)b] 

i<j 

bj = exp(7Ti c\ eje] (12) 

i<j 


so that we obtain a fermionic chain model. Since only 
nearest neighbor hopping is allowed the Hamiltonian re¬ 
mains bilinear when expressed in fermions. The grand 
canonical partition function of these fermions is simple 
and given by 


z g{P,Ha) = Tr [exp (—/3F 0 +/3 a*o™)] (13) 

where fi = . fij counts the number of fermions and 

(n) = N is the expectation value for the number of ther¬ 
mally excited fermions in the chain. This is the total 
number of thermally excited bosons irrespective of their 
flavor. The Hamiltonian Hq in terms of fermions is given 

by 

H 0 = Y2 c i c k uj ( k )’ ( i4 ) 

k 


where the Fourier transformation to momentum space 
has already been carried out to diagonalize the Hamilto¬ 
nian. The grand canonical partition function in the eigen 
basis of the fermions reads 

Zg{P, Ho) = I| + ex P (“Z 3 “ Mo])) (15) 

k 


implying the grand canonical potential 


H = Y2, f A + Y^ &t a b i a MP,Ho) = -^-\nZ g (P,po) (16a) 

7=1, a x 

i=L w = — Vln[l +exp(-/3w(fc))exp(^/r 0 )]. 

-Ei(iw+ h 4 (io) 0 ■ 

i=l,a 


(16b) 
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In the grand canonical ensemble, the particle number N 
is no canonic thermodynamic variable so that we cannot 
simply add the effect of different flavors. Hence we per¬ 
form a Legendre transformation to the free energy 

F o (0,N) = J + n 0 N (17a) 

dF 0 = —SdT + n 0 d N. (17b) 


which reads 

Fo(P,N) = —!-^Jln[l + exp(-^w(fc))exp(/3/x 0 (AT))] 

k 

+ Fo (N)N. (17c) 

At fixed N the effect of the different flavors can be ac¬ 
counted for. The additional degree of freedom leads to 
an entropy increase and hence to a decrease in the free 
energy 


F = F 0 - T Sfl avor (18a) 

= F 0 - i In (N?) (18b) 

where the flavor entropy S'flavor = In is determined 

by the number of flavors Nf and the number of thermal 
excitations N. The free energy F takes the effect of dif¬ 
ferent flavors into account. 

In a next step, we calculate the modified chemical poten¬ 
tial 


fi = d N F(p,N) (19a) 

= no - ^ In (N { ). (19b) 

vf(P) 

The effect of different flavors is seen in the increase of the 
effective chemical potential. A short calculation provides 
us with the modified grand canonical potential J = F — 
fiN from which we finally obtain the thermal occupation 


nexact(T) = ^ (20a) 

= i^[exp(/lw(fc))+ 7Vf] _1 . (20b) 

k 

This result shows the correct behaviour for T = 0 =>■ 
n(T) = 0 and for T —> oo => n(T) = l/(l+lVf). For Nf = 
1 the Fermi-function is recovered as it has to be. 

We stress that the mapping to free fermions by the 
Jordan-Wigner transformation leads to a separation of 
the flavor and the particle dynamics only for a fixed num¬ 
ber of particles. If we consider dynamic correlation func¬ 
tions which involve changes of the number of bosons we 
do not know of a mapping to a solvable model which 
makes these correlations tractable. But further attempts 
are called for. 


III. Diagrammatic Approach 

In this section, we introduce the diagrammatic ap¬ 
proach to treat the hardcore repulsion at finite temper¬ 
atures. The aim is an expansion in the small param¬ 
eter exp(— 0A). The approach has been already intro¬ 
duced in [44J for bosons of a single flavor. Here we focus 
on the changes necessary to treat multi-flavor hardcore 
bosons. 

There is no diagrammatic perturbation theory for hard¬ 
core bosons at finite temperatures. Thus, the key idea 
is to start from normal bosons and to enforce the hard¬ 
core constraint by means of an infinite on-site interaction 
U —> oo. For nuclear matter and He 3 this idea was in¬ 
troduced under the name of Bruckner theory [45]. It is 
important to see the difference of our approach to the 
one in Refs. [45j] and (47l-[49j where the leading effects of 
quantum fluctuations are addressed by Bruckner theory. 
In the present application, the Bruckner theory captures 
the leading thermal effects on the spectrum, e.g., the 
thermal broadening of single-particle lines. 

We replace the initial Hamiltonian FI o in Eq. © by 


H = H 0 + Hu (21a) 

Hu = jim b lAAA«’ ( 21b ) 

i cn,(p 

where U is the strength of the auxiliary local repulsion. 
The main difference to [44| is that different flavors repel 
each other on the same site as well. 

We want to calculate the single particle spectral function 
A(p,u>). At zero temperature, the spectral function is a 
(5-function located at the energy given by the dispersion, 
i.e., A(p,u)) = S(co — oj(p))- At finite temperature the 
excitations are scattered by the thermally excited back¬ 
ground generating a broadening of the spectral function 
in frequency space. 

Since we have to take the limit U —> oo, truncated per¬ 
turbation theory in the repulsion strength U is not use¬ 
ful. Instead, we aim at a low-temperature expansion. 
On the diagrammatic level one observes that each prop¬ 
agator running backwards in imaginary time implies a 
factor exp(—/3A). Thus, in first order in exp(—/3A) one 
has to sum all self-energy diagrams with one backward 
running propagator. This amounts up to the summation 
of the ladder diagrams as depicted in Figs. |T] and [2j The 
gray box indicates the scattering amplitude or effective 
interaction T given in Fig. [2] 

The scattering amplitude fulfills the Bethe-Salpeter equa¬ 
tion given in Fig. [3] which is essentially a geometric series. 
Due to the simple structure of the local interaction, it can 
be easily evaluated leading to 


r «A(p) 


u 

_ L£ _ 

1 + TbE G a (P + Q)G*{-Q)' 
Q 


( 22 ) 
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Figure 1. Self-energy diagrams in leading order in exp(— {3A). 
The first diagram are Hartree-like diagrams and the second 
refers to the Fock-like diagrams with a renormalized interac¬ 
tion given by the scattering amplitude T. 



Figure 2. Definition of the scattering amplitude T 


The capital letters like P and Q denote the 2-moment 
P = ( p , iojp ) and so on. In contrast to the situation in our 
previous work 0 ] the scattering amplitude depends on 
two flavor indices a and (f> denoting the different flavors 
of the propagators G in Eq. (El) . The same applies to 
the spectral functions p p (x) and p p (x) introduced in Ref. 
0 . 

P + K P + Q P + K P + Q' P + Q 


+ £ 

-K -Q -K -Q' -Q 

Figure 3. Bethe-Salpeter equation for the scattering ampli¬ 
tude 

Evaluating the diagrams in Fig. [T] yields the expres¬ 
sion 

s “< p ) = z + a <+(P + K), (23) 

0 K 

for the self-energy. The term proportional to 1 stems 
from the Hartree-like diagram whereas the term pro¬ 
portional to Sarf stems from the Fock-like diagram. If 
the flavors are degenerate, i.e., the propagator G a (P) = 
G^(P) = G(P) does not depend on the index a , one can 
simply replace the sum over the flavor index 4> in (12,311 by 
a factor 

£“(P) = 1 + N{ ^ G a (K)T a ’ a (P + K). (24) 

K 

Once the self-energy is calculated, one determines the 
spectral function 

A a (p,u>) = 

— 1 Im£ a (oc,p) 

71 (oj — ui(p) — Re£“(w,p )) 2 + (Im£ a (w,p )) 2 

(25) 

We conclude that the main difference to the case of hard¬ 
core bosons of a single flavor is the factor 1 + Ay in Eq. 
(El) for degenerate flavors. 


We emphasize that the spectral function is computed self- 
consistently, i.e., we use the bare propagators Gq(P) = 
1 /(iuj p —uj a (p)) only as initial guess for the dressed Green 
function G a (P ) in the diagrams. After the first compu¬ 
tation of the self-energy we iterate all following calcula¬ 
tions with the full propagator G a (P) = l/(icu p — uj a (p) — 
£“(P)) with the previously obtained self-energy until we 
reach self-consistency within numerical tolerance. 

For normal bosons, i.e., without infinite interactions, the 
self-consistent approximation based on ladder diagrams 
is a conser ving ap proximation as introduced by Kadanoff 
and Baym [53|,|54|. Such approximations are derived from 
a generating functional $ depending on the dressed prop¬ 
agators. The self-energy diagrams are obtained from <f> 
by functional derivation with respect to the propagators. 
This construction ensure that conservation laws and sum 
rules are fulfilled. Thus the question arises whether the 
approach is also conserving for hardcore bosons. 
Kadanoff and Baym showed the applicability of the T- 
matrix approximation to single-flavor bosons with hard¬ 
core interactions f46j in continuum models. They also 
established that the T-matrix approximation is conserv¬ 
ing. To our knowledge, however, they did not establish 
generally that the T-matrix approximation is conserving 
for infinite interactions. In Ref. 0 there is no discussion 
of the difference between first integrating over frequency 
and then taking the limit U —> oo and the opposite se¬ 
quence. 

Empirically, we found previously j0 that our approach 
conserves the sum rules for hardcore bosons of a single fla¬ 
vor within numerical accuracy. We emphasize that this is 
a non-trivial statement which does not result from con¬ 
tinuity in the limit U —> oo. The integration over all 
frequencies from — oo to oo does not commute with the 
limit U —> oo because the latter shifts weight to infinity. 
This weight is lost in the integration explaining why the 
weight of the spectral functions of hardcore bosons is re¬ 
duced as expressed by the right hand side of Eq. ©. The 
other sequence, namely integrating first and then having 
U tend to infinity describes normal bosons and does not 
alter the sum rule which remains unity. We will show 
this explicitly in Section fVl 

IV. Results for Multi-Flavor Bosons 

In this section we discuss the results of the Bruckner ap¬ 
proach in the multi-flavor case. We focus on the changes 
upon passing from Nf = 1 to Nf = 2 and Nf = 3. Es¬ 
pecially the case TVf = 3 is relevant because it describes 
triplon excitations in dimerized spin systems. We will 
evaluate sum rules and the thermal occupation in detail 
to assess the validity of the Briickner approach on the 
level of ladder diagrams at finite temperature. 

All results are obtained for the nearest neighbor model 
in Eq. COD without any additional interactions besides 
the hardcore constraint. We primarily show results for 
W/A = 0.5, i.e., for a narrow band. For comparison, 
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Figure 4. Thermal occupation n(T) in the narrow-band case 
W = 0.5A for various numbers of flavors Nf as function of 
temperature. Comparison between the exact curves (solid 
lines) and the diagrammatic results (dashed lines). 



T[A] 


Figure 5. Thermal occupation n(T) in the wide-band case 
W = 4A for various numbers of flavors Nf as function of 
temperature. Comparison between the exact curves (solid 
lines) and the diagrammatic results (dashed lines). 

we will also show the thermal occupation function in a 
wide-band case, W /A = 4. 

First, we study the thermal occupation which can be 
computed from the spectral function by 

2-7T 27T OO 

n(r) = h J= JJ Y dujdk ■ 

0 0 —oo 

(26) 

Note that the thermal occupation is independent of the 
index a because the bosons of different flavor are de¬ 
generate. Figs. [I] and [5] depict the results for different 
numbers of flavors obtained by the Bruckner approach 
and compare them to the exact expression (I20bl) . 

Since our approach is based on a low-temperature ap¬ 


proximation in the expansion parameter exp(— /3A) we 
expect deviations to occur if the temperature approaches 
the value of the energy gap A. Interestingly, the Nf = 1 
thermal occupation is approximated very well even at el¬ 
evated temperatures 0] and even for vanishing gap im, 
only the relative error diverges for A —>■ 0. This holds for 
n(T); the deviations in the s pec tral function A(p, w) are 
larger, see Section V in Ref. |44| . 

As expected, the diagrammatic and the exact results in 
Figs. Q] and [5] agree exellcntly for small temperatures T < 
0.3A. For large temperatures T —> oo the exact curves 
yield the value n(T oo) = 1/(1 + Nf) corresponding to 
the value in isolated dimers where all states are equally 
probable. In contrast to the Nf = 1 case, discernible 
deviations occur between the approximate and the exact 
curves for Nf > 1 upon increasing temperature. The 
Bruckner approach on the level of ladder diagrams for 
Nf > 1 tends to underestimate the thermal occupation 
at higher temperatures. But we emphasize that for low 
temperatures the approach works as expected because it 
is exact in order exp(—/3A), see below. 

There is also a significant difference between the perfor¬ 
mance of the ladder approximation for a narrow and for 
a wide band. At given gap and temperature the relative 
deviations are significantly smaller for wider bands. In 
the W = 4A case, the diagrammatic curves stay close 
to the exact curves even for T > 0.6A. In stark con¬ 
trast, for W = 0.5A, Nf = 3, and T > 0.75A the ap¬ 
proximate thermal occupation even shows an unphysical 
non-monotonic behavior as function of T. 

The difference between the two cases can be explained by 
the different total occupations at the same temperature. 
For a wide band the fraction of the Brillouin zone with 
low-lying excitations in the range of T is smaller than 
for a narrow band so that the total number of thermally 
excited particles is significantly lower. For this reason, 
the corrections to the scattering processes considered by 
the ladder diagrams are less important. 

This observation implies an important message also for 
the applicability of the ladder approximation in higher 
dimensions. There, the fraction of the Brillouin zone at 
low energies for a given gap is smaller than in one di¬ 
mension: Assume the gap to occur at k = 0 then all 
components of k must be small for the total energy to 
be small. For the same parameters temperature T, gap 
A, and band width W one may roughly estimate that 
the thermal occupation in d dimensions nd(T) scales like 
n\(T) d . Thus we expect that for given parameters the 
accuracy of the approximation on the level of ladder dia¬ 
grams is significantly higher in higher dimensions. 

Next, we will investigate the deviations in one dimension 
in more detail. To this end, we study the sum rule of 
the spectral function A(p,oj). In Eq. 121) we saw already 
that the integration of A(p,u>) over frequency ui is con¬ 
nected to the thermal occupation. Therefore we define 




















Figure 6. The sum rule R[n(T)\ from Eq. (1271) for vari¬ 
ous numbers of flavors Nf as function of temperature. The 
bandwidth is fixed to W = 0.5A. We compare the quantity 
R[n(T)] using the diagrammatic n(T) obtained from Eq. (1261) 
and f?[n e xact (T)] using the exact result from Eq. (120bb . 

the expression 

OO 

R[n(T)]:= J A a {x)dx + (1 + N t )n(T), (27) 

— OO 

which should be equal to unity for the exact spec¬ 
tral function. For brevity, we have introduced A(x) = 
j- A a (p, x), i.e., the spectral function averaged over 

p 

momentum p. 

The sum rule © states that the integral over the momen¬ 
tum dependent spectral function A(p, cu) is momentum 
independent. Any approximation to calculate the spec¬ 
tral function is prone to spoil this feature. But we observe 
in our numerical data that the Briickner approach com¬ 
pletely conserves the exact p-independence within nu¬ 
merical accuracy for all numbers of flavors Nf. Thus, 
there is no need to discuss a momentum dependence of 
the sum rule. It is sufficient to analyse the quantity 
R[n(T)}. 

We plot R[n(T)] as function of temperature in Fig. [G] 
Since for Nf > 1 the thermal occupation n(T) obtained 
from the spectral function differs from the exact expres¬ 
sion in Eq. (I20bl) we depict results for (1271) for both occu¬ 
pation functions, the one determined diagrammatically 
and the exact one. 

For Nf = 1 the sum rule is fulfilled within numerical ac¬ 
curacy even for temperatures close to the gap A. This 
changes for higher flavor numbers Nf > 1. At low tem¬ 
peratures, the sum rule is well fulfilled because by con¬ 
struction the approach is exact in order exp(—/3A). But 
deviations become discernible at higher temperatures. 
Some spectral weight seems to be missing in the approx¬ 
imate Bruckner approach. Note that the sum rule turns 
out to be better fulfilled using the exact n exa , c t{T) than 
using the diagrammatically obtained n(T). Comparing 



Figure 7. Quantitative analysis of the deviation in the calcula¬ 
tions of the sum rule for Nf > 1 in Fig. [S] as function of inverse 
temperature /3A. The bandwidth is fixed to W = 0.5A. The 
y -axis shows the deviation of the sum rule from 1 using the di¬ 
agrammatic thermal occupation n(T) from Eq. (1261) . Straight 
lines indicate exponential dependence of the inverse temper¬ 
ature as expected. The slope indicates the prefactor in the 
argument of the exponential. 

the case Nf = 2 and Nf = 3, one realizes that the rel¬ 
ative deviation increases upon increasing the number of 
flavors. 

For a quantitative analysis of the deviations for Nf > 1, 
Fig. 0 depicts the logarithmic deviation In 11 — i?[n(T)]| 
of the sum rule in Fig. [5] as function of inverse temper¬ 
ature /3A. We use the diagrammatically obtained occu¬ 
pation function n(T) from (1261) . For low temperatures 
T < 0.3A, i.e., high inverse temperature /3A > 3.3, the 
deviations smaller than the numerical accuracy of the 
data 

Upon increasing temperature the logarithmic deviation 
scales linearly with inverse temperature indicating an ex¬ 
ponential dependence on inverse temperature. We fit¬ 
ted the functions Aexp(—/3A) and Bex p(—2/3A) to the 
data in the low temperature regime fd A > 3. From 
these fits it is evident that the deviation indeed scales 
with exp(—2/3A). This in accord with the construction 
of the approximation which was designed to be exact 
in linear order in exp(—/3A). Thus in the next order, 
i.e., oc exp(—2/3A), deviations cannot be excluded. The 
analysis provided underpins the correct reasoning in the 
choice of diagrams and in the implementation of the nu¬ 
merical calculations. 

In Fig. 0 we additionally show the logarithmic deviation 
of the sum rule using the exact occupation function n(T) 
from Eq. (I20bl) . Again the deviation scales as exp(—2/3A) 
indicating that the integral over the spectral function is 
only correct in first order. 

Next, we compare the thermal occupation obtained from 
the ladder approximation in Eq. m to the exact expres¬ 
sion from Eq. (120bl) in Fig. [9] Again we show the logarith¬ 
mic difference of the two expressions as function of the 
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Figure 8. Quantitative analysis of the deviation in the calcula¬ 
tions of the sum rule for Nf > 1 in Fig. HU as function of inverse 
temperature /3A. The bandwidth is fixed to W = 0.5A. The 
y- axis shows the deviation of the sum rule from 1 using the ex¬ 
act thermal occupation n exac t(T) from Eq. (I20bll . Otherwise 
as Fig. [7] 

inverse temperature. Surprisingly, the difference scales 
like exp(—3/3A). This indicates a result better than one 
could expect naively from the construction of the ladder 
approximation. 

We attribute this better-than-expected agreement to 
the additional averaging of the spectral function over 
the Bose function in Eq. (1261) . This appears to cor¬ 
rect the second order deviation exp(—2/3A). We refer 
the reader to the analogy in leading order: The zero- 
temperature spectral function corresponds to the zeroth 
order exp(—0/JA) and leads via the Bose function to the 
correct result for the thermal occupation in first order 
exp(— pA). Thus one can expect that the spectral func¬ 
tion in first order implies the correct thermal occupation 
in second order. 


V. Analysis at Finite U 

In the previous section, we saw that for Nf = 1 the 
sum rule is conserved within numerical accuracy and 
the thermal occupation is approximated very well by 
the Bruckner approach based on ladder diagrams. In 
contrast, the multi-flavor case represents a less power¬ 
ful approximation. But we verified that the occurring 
deviations are consistent with the systematic construc¬ 
tion of the approximation as an expansion in powers of 
exp(— /3A). At low temperatures, the approach works as 
expected. 

Still, these findings are puzzling in view of the concept 
of conser ving ap proximations introduced by Baym and 
Kadanoff [53l l54l | for normal bosons and fermions. If the 
choice of diagrams yields a conserving approximation for 
Nf = 1 it is unexpected that this is no longer the case for 
Nf > 1 . We stated, however, already above that the case 


Figure 9. Quantitative analysis of the deviations in the cal¬ 
culations of the thermal occupation in Fig. [5] as function of 
inverse temperature /3A for Nf > 1 . The bandwidth is fixed 
to W = 0.5A. The j/-axis displays the difference between 
Uexact(T) from Eq. (I20bl) and the diagrammatic n(T ) from 
Eq. (HU. Unexpectedly, one finds that the deviations scale 
proportional to exp(— 3/3A). 

of hardcore bosons is more subtle anyway. To our knowl¬ 
edge, there is no general argument to prove that a set of 
diagrams yielding a conserving approximation for normal 
bosons also implies a conserving approximation for hard¬ 
core bosons due to the non-trivial sequence of integration 
over all frequencies and the limit U —> oo. 

The fundamental question is whether the concept of 
conserving approximations is generally transferable from 
normal bosons to hardcore bosons, i.e., independent of 
the flavor number Nf. To understand this issue better 
and to verify that our choice of diagrams provides a con¬ 
serving approximation for normal bosons we analyze the 
approximation of ladder diagrams for finite repulsion U. 
For finite repulsion, the particles considered are normal 
strongly interacting bosons. Here the Baym-Kadanoff 
concept of conserving approximations is applicable. We 
study how the sum rules behave in this case and upon 
the limit U —> oo. 

In Figs. [TO] and [TT] we compare the spectral line of the 
scattered boson for U = oo with results for U — 5A, 
U = 10A and U = 20A . For low temperatures, already 
the case U = 5A approximates the U = oo result quite 
well. Upon increasing temperature larger deviations oc¬ 
cur. But clearly, the data for finite values of U slowly 
converge towards to the results at U = oo. 

For bosons with finite on-site interaction the general sum 
rule in Eq. ([8]) yields unity for all temperatures 

OO 

J ^Bosonic (p,w)du; = 1. (28) 

— OO 

We carefully checked that the approximation based on 
ladder diagrams indeed fulfills this sum rule independent 
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Figure 10. Comparison of the spectral lines A(p, to) at the 
gap mode p = 0 and at the maximum mode p = n for various 
values of the on-site repulsion U. The temperature is fixed to 
T = 0.3A and the number of flavors is Nf = 3. 


Figure 12. Full frequency range comparison of the spectral 
function A(p,u)) of the gap mode p = 0 for various values of 
the on-site repulsion U. The temperature is fixed to T = 0.5A 
and the number of flavors is Nf = 3. At negative frequencies 
there is only negligible weight for the considered parameters. 
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Figure 11. Comparison of the spectral lines A(p, u>) at the 
gap mode p = 0 and at the maximum mode p = n for various 
values of the on-site repulsion U. The temperature is fixed to 
T = 0.6A and the number of flavors is Nf = 3. 


of the number of flavors Nf and of the temperatures T 
within numerical accuracy. Since the spectral lines at 
finite U converge to the result at U = oo on increasing 
U, there must be additional weight somewhere in the 
spectral function to explain the difference to the spectral 
weight in the U = oo case. 

We find the additional weight at energies in the range 
w « U. Fig. [l2] displays the spectral function A(p, oj) 
over the whole energy range. For finite U, the addi¬ 
tional weight at high energies can be explained by the 
formation of an anti-bound state of two bosons which re¬ 
pel each other on the same site. The signature of this 
anti-bound state is found in the spectral function of the 
scattering amplitude resulting from the Bethe-Salpeter 
equation (not shown). Clearly, the Bethe-Salpeter equa¬ 
tion describes the repeated scattering of two propagating 


particles. The energy of the anti-bound state is of the or¬ 
der of the repulsion U. 

The additional convolution with the third propagator 
only broadens the line by about the band width W. This 
is what is seen in Fig. [T2] at high energies. Note that the 
position of this high-energy feature scales with U. It is 
this feature which contains the missing spectral weight 
for the total sum rule of the single-particle spectral func¬ 
tion of normal bosons. 

The first key observation for general Nf in the above 
data is the confirmation that the two limits U —> oo 
of the Bruckner approach and w —> oo in the integration 
over the frequencies do not commute. The properties of 
the hardcore bosons are retrieved if the limit of infinite 
repulsion is performed first. 

For the finite U case, the ladder diagrams represent a con¬ 
serving approximation independent of Nf. In the infinite 
U case, the approximation appears to be a conserving 
approximation only in the single-flavor case Nf = 1, but 
not for Nf > 1. 

We stress that this does not contradict the results by 
Baym and Kadanoff, because they formulated conserv¬ 
ing approximations for normal bosonic systems with ar¬ 
bitrary, but finite interactions [46j. They also discussed 
that infinite interactions can be treated by the T-matrix 
approximation, but did not consider the subtleties of the 
sequence of limits, i.e., they did not elucidate whether 
this approximation is conserving for infinite interactions. 
The single flavor case suggested that the Baym-Kadanoff 
construction also yields conserving approximations for 
hardcore bosons. But our calculations for the Nf > 1 
falsify this hypothesis. From this observation we deduce 
our second key conclusion that the concept of conserving 
approximations in the sense of Baym and Kadanoff does 
not carry over to hardcore bosons in general. 

The fundamental questions remains why hardcore bosons 
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with multiple flavors behave differently from single-flavor 
hardcore bosons in one dimension. Answering this ques¬ 
tion is beyond the scope of this article although we for¬ 
mulate a hypothesis in the Conclusions [V]] but further 
studies are called for. 


VI. Conclusions 

In this paper, we extended the Bruckner approach from 
the single-flavor case to the multi-flavor case of hardcore 
bosons. Such multi-flavored excitations are for instance 
experimentally relevant in coupled spin dimers systems 
where the excitations are threefold degenerate triplons. 
We found that the required modifications in the ladder 
diagrams for degenerate flavors are minimal. Only the 
prefactor in the self-energy is modified. 

The focus of the present work is to gauge the diagram¬ 
matic approach in the multi-flavor case against exact re¬ 
sults and to study conceptual issues concerning conserv¬ 
ing approximations for hardcore bosons. 

Rigorous results were derived for this purpose for the 
thermal occupation of hardcore bosons in the case of 
a one-dimensional chain with nearest-neighbor hopping. 
The established rigorous relation is of interest in itself 
because it can be used as benchmark for many approxi¬ 
mate analytic or numeric treatments. It could be derived 
by observing that the particle degree of freedom and the 
flavor degree of freedom represent two independent sub¬ 
systems for fixed number of bosons. This helps to treat 
the thermodynamics. But the dynamic correlations, for 
instance the hardcore boson propagator, cannot be dealt 
with because they involve changes of the particle num¬ 
ber. 

Our first key results consists in the benchmarks which 
clearly showed that the multi-flavor approach based on 
ladder diagrams works very well for low temperatures. 
This is expected from the construction of the approach 
as the leading linear order of an expansion in powers of 
exp(—/3A). Thus it can be used up to T « A/2 for nar¬ 
row bands and up to T « A for wide bands. For bosons 
of a single flavor it worked even better as established pre¬ 
viously [13|. At present, we can only speculate why the 
single flavor case is more easily tractable. 

The conceptual issue of conserving approximations con¬ 
sists in the question whether the construction of con¬ 
serving approximations bv Baym and Kadanoff for nor¬ 
mal bosons and fermions [53, [54| carries over to hardcore 
bosons. In the present article, we confirmed and veri¬ 
fied for finite U that our choice of self-consistent ladder 
diagrams represents a conserving approximation in the 
Baym-Kadanoff sense. This holds true independent of 
the flavor number Nf. 

Our previous empirical finding (44| in the single-flavor 
case Nf = 1, indicated that the hardcore boson sum rules 
are fulfilled as well within the numerical accuracy. Thus 


one could speculate that the Baym-Kadanoff construc¬ 
tion carries over to hardcore bosons. But the results of 
the present article for multi-flavor bosons Nf > 1 un¬ 
equivocally show that this is generally not the case. This 
constitutes the second key result of our study. For sev¬ 
eral flavors the sum rules are not conserved by the ladder 
diagram. We stress, however, that the deviations only 
become sizable at higher temperatures. 

We presume that the qualitative difference Nf = 1 and 
Nf > 1 is due to the fact that the case of a single hardcore 
boson in one dimension can be exactly mapped to nor¬ 
mal fermions for which conserving approximations exist. 
In the case of multiple flavors the corresponding Jordan- 
Wigner mapping to multi-flavor fermions does not lead 
to normal fermions because the ensuing fermions are still 
hardcore particles repelling each other infinitely strongly 
if they are of different flavor. We can still use this map¬ 
ping to describe the thermodynamics of multi-flavored 
hardcore bosons for a fixed number of particles. But dy¬ 
namic correlations changing the number of particles in 
the system are not trivially accessible in this way. 

We summarize that the Bruckner approach based on self- 
consistent ladder diagrams also works for multi-flavor 
excitations in the interesting low-temperature regime 
though it is less powerful in comparison to the single- 
flavor case. In particular for wider bands, the approach 
works very well and we argued that the approximation 
should perform even better in higher dimensions, where 
the diagrammatic approach can be used as well without 
conceptual obstacle. 

The presented analysis paves the way for the treatment 
of more realistic models, i.e., models closer to experimen¬ 
tally relevant systems. Obviously, the dispersion can be 
straightforwardly generalized. But also further interac¬ 
tions beyond the hardcore repulsion need to be included. 
The analysis of higher dimensional models is called for 
and first results in the single-flavor case have become 
available very recently [sTJ] . 

Having these extensions in mind, we stress that the 
Bruckner approach is a promising tool to describe exper¬ 
imental observations, such as asymmetric thermal line- 
shape broadening as well as thermal narrowing of the 
single excitation bandwidth. 

Further conceptual extensions concern the addition of the 
considered diagrams to those with two backward running 
propagators providing an approximation exact up to or¬ 
der exp(— 2/3 A). 
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